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ABSTRACT 

After a careful analysis of the instrumental effects on the Poisson noise to demonstrate the feasibility 
of detailed stochastic variability studies with the Swift X-Ray Telescope (XRT), we analyze the vari- 
ability of the black hole X-ray binary SWIFT J1753. 5-0127 in all XRT observations during 2005-2010. 
We present the evolution of the power spectral components along the outburst in two energy bands: 
soft (0.5-2 keV) and hard (2-10 keV), and in the hard band find results consistent with those from 
the Rossi X-ray Timing Explorer (RXTE). The advantage of the XRT is that we can also explore the 
soft band not covered by RXTE. The source has previously been suggested to host an accretion disk 
extending down to close to the black hole in the low hard state, and to show low frequency variability 
in the soft band intrinsic to this disk. Our results are consistent with this, with at low intensities 
stronger low-frequency variability in the soft than in the hard band. From our analysis we are able to 
present the first measurements of the soft band variability in the peak of the outburst and find it to be 
less variable than the hard band, especially at high frequencies, opposite to what is seen at low inten- 
sity. Both results can be explained within the framework of a simple two emission-region model where 
the hot flow is more variable in the peak of the outburst and the disk is more variable at low intensities. 

Subject headings: X-ray binaries: general — black hole candidate : individual SWIFT J1753. 5-0127 



1. INTRODUCTION 

Transient stellar-mass black hole candidate X-ray 
binaries (BHBs) in outburst broadly exhibit two dif- 
ferent states: low luminosity hard states in which the 
spectrum is dominated by hard power law emission 
(out to a few tens of keV) from a hot inner flow, inner 
corona, and/or the base of a jet, and high luminosity 
soft states in which the spectrum is dominated by 
black-body emission (peaking at up to a few keV) from 
the accretion disk, often accompanied by a power law 
tail. There is no agreement about the exact structure 
and origin of the power law emitting region, nor on the 
disk geometry: while previously it was thought that in 
the hard state the disk is truncated at large (at least 
several 10i? g ) radii, recently it has been debated if 
instead it extends all the way down to the innermost 
stable circular orbit (ISCO). See Belloni et al. (2005), 
van der Klis (2006) and Remillard & McClintock (2006) 
for detailed descriptions of the different states along an 
outburst and Done et al. (2007) for a detailed discussion 
of the different models. 

The ubiquitous X-ray variability in BHBs could 
provide useful constraints on the structure and behavior 
of the emitting regions, beyond that provided by the 
X-ray spectrum. Lyubarskii (1997) proposed a model 
which is becoming widely accepted (e.g. Churazov 
et al. 2001; Uttley et al. 2005; Done et al. 2007), where 
the variability originates as mass-accretion fluctuations 
in the accretion flow, which propagate towards and 
through the X-ray emitting regions to produce the 
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observed flux variability. Since the variability in the 
hard states is strong (few tens of percent fractional 
rms) while in the soft states it is much weaker (a few 
percent), Churazov et al. (2001) suggested that the hot, 
power law emitting inner flow is responsible for this 
variability, while the cool disk is stable. Furthermore, 
they showed that in a 'standard' geometrically-thin 
cool disk, high frequency fluctuations can survive only 
if they arise at small radii, as at larger radii, viscous 
damping will prevent them from propagating inward. 
In the hard state, the power spectrum consists of broad 
band-limited noise accompanied by Quasi Periodic 
Oscillations (QPOs). The low frequency break in the 
noise power spectrum could be associated with the 
transition radius to the hot flow at the inner edge of 
the cool disk, the QPO could be the Lense-Thirring 
precession of the hot flow, and the high frequency 
break may be associated with the innermost radius of 
the hot flow (Stella & Vietri 1998, Churazov et al. 2001). 

The Rossi X-ray Timing Explorer (RXTE) Pro- 
portional Counter Array (PCA) was sensitive in the > 2 
keV band dominated by the power law emission. Most 
variability studies until recently were done with RXTE, 
so the evolution of the variability below 2 keV remained 
unexplored. Recent work by Wilkinson & Uttley (2009) 
and Uttley et al. (2011) using XMM-Newton showed 
that in the hard state of the BHBs SWIFT J1753.5-0127, 
GX 339-4 and Cyg X-l the disk contributes significantly 
to variability in the soft band (< 2 keV) at Fourier 
frequencies below 1 Hz and moreover, the disk variations 
precede the correlated power law variations by a few 
tenths of a second. These results cast doubt on the 
notion of a stable disk not contributing to the source 
variability, implying that (at least in the hard state and 
at low frequencies) substantial flux variability is con- 
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Fig. 1. — Light curve of the outburst during 2005-2010 (top) and 
corresponding hardness-intensity diagram (bottom) in the energy 
bands indicated. Each point represents average intensity and hard 
color per observation. The standard deviations are up to 8 c/s for 
the intensity and 0.05 for the hard color. No corrections for pilc-up 
or background were made. The horizontal lines at 150 c/s indicate 
the level above which pile-up affects the data. 

tributed by the cool disk, highlighting the importance 
of exploring the variability in the soft band. Like XMM- 
Newton, the Swift X-ray Telescope (XRT) has a lower 
energy bound (0.3 keV) than the RXTE PCA, providing 
another opportunity to explore the variability at lower 
energies dominated by the disk emission, and to di- 
rectly address the question of the origin of the variability. 

SWIFT J1753. 5-0127 was discovered with the Swift 
Burst Alert Telescope on 2005 May 30 (Palmer et al. 
2005) and has been in outburst since then. The source 
remained in the hard state (Miller et al. 2006 and 
Cadolle Bel et al. 2007, Zhang et al. 2007, Ramadevi 
& Seetha 2007, Chiang et al. 2010, henceforth C07, 
Z07, RS07 and ChlO, respectively) until June 2009 
(MJD 55010) but then was in a softer state until at 
least July 2010 (MJD 55404) (Soleri et al. 2012). 
During 2005-06, RXTE observations show white noise 
(frequency independent) at low (0.2-0.4 Hz) frequencies, 
a QPO (~0.7 Hz) and red noise (power decreases with 
frequency) at high (1.5—2 Hz) frequencies, which is 
typical of the hard state (Z07, RS07). Spectral analyses 
with XMM-Newton and RXTE suggested the presence 
of a cool accretion disk extending down to near the 
ISCO that could be modelled with a disk black-body 



component in addition to the power law (Miller et al. 
2006, ChlO). This is at variance with the picture in 
which the disk is truncated at larger radii in the hard 
state. There are also arguments against the presence of 
a soft disk- like component (see, e.g., Done et al. 2007 
and Hiemstra et al. 2009), which, however, do not 
refute the evidence for intrinsic disk variability found by 
Wilkinson & Uttley (2009). 

In this work, we use the Swift XRT to study the 
energy-dependent behavior of the variability of the BHB 
SWIFT J1753. 5-0127. While some XRT power spectra 
of BHBs were previously reported (Kennea et al. 2011, 
Curran et al. 2011), the present work is the first study 
of the energy dependence extending below 2 keV of all 
power spectral components in a BHB. As a prerequisite 
to this work we study the instrumental effects on the 
Poisson-noise spectrum (see Appendix) and demonstrate 
that detailed variability studies are feasible with Swift 
XRT. In Section 2, we describe the observations and data 
analysis. In Section 3.1, we discuss the evolution of the 
outburst using light curves and the hardness intensity 
diagram. In Section 3.2, we first investigate the energy 
dependence of the variability in a model-independent 
way, and then fit Lorentzians to the power spectra of all 
observations to study the evolution of the components 
and their energy dependence. We also study the energy 
dependence of the power spectra inside the 0.5-2 keV 
band. Finally, in Section 3.3 we study the correlations 
between the variability parameters and source intensity. 
In Section 4, we discuss the interpretation of our results, 
and their model implications. 

2. OBSERVATIONS AND DATA ANALYSIS 

We analysed all 66 observations taken in the Window 
Timing (WT) mode with the X-Ray Telescope (XRT; 
Burrows et al. 2005) on board the Swift satellite between 
July 1, 2005 and May 21, 2010. Observations lasted 
between 0.1 and 11.0 ksec, and contained between 1 and 
34 Good Time Intervals (GTIs) of 0.1-2.0 ksec. The 
CCD was operated in the WT mode in wt2 configura- 
tion, in which every 1.766 ms a one-dimensional image 
of 200 pixels is read out from the central columns of the 
CCD (see also the Appendix). 

We processed the raw data with the task xrtpipeline 
v0.12.6 using standard quality cuts and only grade 
events (Evans et al. 2007). The data were extracted in 
two ways. In the first extraction, using xselect v2.4b, 
we generated light curves in 0.5-2 and 2-10 keV energy 
bands by adding up the counts in the entire 200-pixel 
image without any pile-up (two photons recorded as a 
single event, see below) or background corrections. 

For the second extraction, aimed at generating the 
power spectra, we first determined the source region on 
the CCD. As there is no point-spread function for the 
WT mode, we fitted a Lorentzian to obtain the source 
centroid. Reduced x 2 values were 0.5-3 (193 degrees of 
freedom), occasionally much higher due to bad pixels; 
in practice this method always sufficed for our purpose. 
The Lorentzian full width at half maximum (FWHM) is 
typically 5 pixels. A 61-pixel (144") source region was 
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Fig. 2.— Light curves of 2005, 2006-08 and 2010 in energy bands 
indicated. Each point represents average intensity per observation. 
The standard deviations are up to 6 c/s. No corrections for pile-up 
or background were made. Note the different scales of the y-axes. 



selected centered on the centroid pixel, plus, contiguous 
with this and on the side of the image that is larger, 
a background region of the same size (or smaller if the 
edge of the image was reached). Count rates exceeding 
150 c/s are at risk of pile- up (Evans et al. 2007). 
To mitigate this, for each GTI the central pixel and, 
if necessary, additional pairs of pixels symmetrically 
around it are removed until the count rate is below 
150 c/s, basically censoring the piled- up data (Evans 
et al. 2007). This introduces artificial jumps in count 
rate between GTIs; the power spectra are not affected 
by this, as they never straddle GTIs. 

Leahy-normalized (Leahy et al. 1983) fast Fourier- 
transform power spectra were calculated of 28.93-s 
continuous intervals containing 16384 1.766-ms samples, 
giving a Nyquist frequency of 283.126 Hz and 8192 
34.6-mHz frequency bins per power spectrum. The 
power spectra were then averaged per Swift observation. 
We also calculated power spectra averaged over multiple 
observations using longer (57.86-s) continuous intervals. 
No background correction was applied prior to the gen- 



eration of the power spectra. To facilitate comparison 
with RXTE, two energy bands were used: 2-10 keV 
(also covered by RXTE) and 0.5-2 keV (not covered 
by RXTE), and occasionally also narrower bands below 
2 keV. 

Due to instrumental effects (see Appendix) , we analysed 
the power spectra in the frequency range <100 Hz only. 
The power spectrum of each observation was fitted with 
a constant for the Poisson level plus several Lorentzians 
in the "f max " representation (Belloni, Psaltis, & van 
der Klis 2002), with as fit parameters the characteristic 

frequency v max = voy/l + 1/(4Q 2 ), the quality factor 
Q= fo/FWHM, and the integrated power P, where 
vq is the centroid frequency and FWHM is the full 
width at half maximum of the Lorentzian. When Q 
turned out negative, it was fixed to (i.e., we fitted a 
zero-centred Lorentzian) ; this did not significantly affect 
the other parameters. We only report components with 
a single-trial significance P/<Jp_ > 3.0, with op_ the 
negative error on P calculated using A\ 2 = 1. 

3. RESULTS 

3.1. Light curves and color diagram 

Figure 1 shows the light curve and hardness-intensity 
diagram (HID) of the 2005-10 observations. Intensity 
is the count rate in the 0.5-10 keV energy band, while 
hardness (or "hard color") is the 2-10/0.5-1 keV count 
rate ratio. Data with intensity over 150 c/s can be 
affected by pile-up and must be interpreted with cau- 
tion. In Figure 2, light curves are presented separately 
for 2005, 2006-08 and 2010 in the 0.5-2 ('soft') and 
2-10 keV ('hard') bands. 

The 2005 light curve was fast-rise exponential-decay, 
consistent with the RXTE results (C07, Z07, RS07 and 
ChlO). The flux peaked at 138 c/s and 50 c/s in the 
soft and hard band, respectively. By MJD 53600 it had 
decayed to below 20 c/s per band, where it remained 
throughout 2006-08. From the peak of the outburst in 
2005 until 2007, the source gradually became harder, 
as can be seen in the HID, with the first observation in 
2007 (MJD 54248.7) the hardest; thereafter it gradually 
softened. When the source was observed again in 2010, 
the intensity had increased to 100 c/s in the soft band 
while in the hard band it was still below 20 c/s. So, 
consistent with the report by Soleri et al., the source was 
much softer, as can also be seen in the HID (hardness < 
1), with the first observation in 2010 (MJD 55313.5) the 
softest. We occasionally refer to the spectral state in 
2010 as soft without implying any particular canonical 
state viz. High Soft State or Very High State. 

3.2. Timing analysis 
3.2.1. Model-independent analysis 

We first investigate the energy dependence of the 
variability in a model-independent way. We calculated 
the fractional rms amplitude from the power spectra 
over the frequency ranges 0.035-0.4 Hz and 0.4-10 Hz 
(chosen based on the frequencies of the power spectral 
components, see below) separately in the hard and 
soft band. Figure 3 shows the rms over these two 
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Fig. 3. — Dependence of fractional rms amplitude on intensity in 
energy bands and frequency ranges as indicated. The points below 
the vertical grey line represent the averages over MJD 53638-53640 
and of 2006, 2007 and 2008; all other points represent individual 
observations in 2005. 

frequency ranges and energy bands as a function of 
intensity for individual observations up to MJD 53638, 
and for the average power spectra of observations 
in MJD 53638-53640 and the years 2006, 2007 and 
2008 (as the count rate was too low to use individual 
observations from MJD 53638 until 2010). The 2010 
data has very little variability (as expected in a soft 
state) which is poorly constrained and hence not plotted. 
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Fig. 4. — The rms normalized power spectrum of an XRT obser- 
vation (00030090015, MJD 53562.339) in the peak of outburst, in 
soft (0.5-2 keV, 91 c/s) and hard (2-10 keV, 38 c/s) bands. The 
best fit model using three Lorentzians is shown in each frame. 
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At low intensity (<40 c/s), the <0.4 Hz variability 
is a bit weaker in the hard band (18-24% rms) than in 
the soft band (20-30%), whereas the >0.4 Hz variability 
has large error bars but is consistent with being equal 
in the two bands at 10-30%. This is consistent with 
what Wilkinson & Uttley (2009) found in the 2006 
XMM-Newton data. 

With Swift we can now report on what happened 
in the peak of the outburst, i.e., at high intensity 
(>40 c/s): the <0.4 Hz variability in the hard band is 
systematically slightly higher than in the soft band for 
intensity >110 c/s (it is similar in both bands between 
40 and 110 c/s). So, going from low to high intensity, 
the behavior swaps between the two energy bands. The 
>0.4 Hz variability above 80 c/s is clearly stronger in 
the hard band (20-25%) than in the soft band (12-18%). 
So, both <0.4 Hz and (using also the information from 
XMM-Newton) >0.4 Hz the fractional rms spectrum 
(as a function of energy) softens as the source decays, 
contrary to the overall source spectrum, which gets 
harder. In terms of differences in fractional rms between 
the peak of the outburst and the decay, the <0.4 Hz 
variability increases as the source decays, but more 
strongly in the soft band than in the hard band. The 
>0.4 Hz variability shows little evidence of change (the 
error bars at low intensity are large), although in the 
soft band there is some evidence for an increase as the 
source drops below 100 c/s. 



Fig. 5. — Intensity, frequency and fractional rms amplitude of 
the components detected in the power spectra in the hard band as 
a function of time in 2005. Each point represents one observation. 
Symbols as indicated. 

3.2.2. Power-spectral fits 

We now investigate how the energy dependence of the 
variability, and in particular the extra variability >0.4 Hz 
we found in the hard band during the peak of the out- 
burst, translates into energy dependencies of the power- 
spectral components. We analyse the power spectra of 
all observations in the hard and soft band. Before we 
study the power spectra in detail, as a sanity check we 
compared our Swift XRT results in the hard band with 
earlier RXTE results. Figure 4, right panel shows a 
representative hard-band power spectrum of an observa- 
tion in the peak of the outburst. Three components are 
detected, which in individual observations during 2005 
vary in characteristic frequency over the ranges 0.07- 
0.3 Hz (at Q=0-2.1), 0.5-0.9 Hz (Q=l-ll) and 0.8- 
1.6 Hz (Q=0-0.8). Henceforth we refer to them as the 
0.1 Hz component, the QPO and the 1 Hz component, 
respectively. Figure 5 shows the light curve, and the 
evolution of frequency and amplitude of the components 
during 2005 in the hard band. The same three compo- 
nents we see with XRT have also been observed with 
RXTE (Z07, RS07, ChlO). As only the frequency evolu- 
tion of the QPO was presented in Z07, RS07 and ChlO, 
we cannot check the other parameters, but our QPO fre- 
quencies match well to the RXTE ones. 
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Fig. 6. — Intensity, frequency and fractional rms amplitude of 
the components detected in the power spectra in soft band as a 
function of time in 2005. Each point represents one observation. 
Symbols as indicated. 
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Fig. 8. — Intensity, frequency and fractional rms amplitude of 
the components detected in the power spectra in the soft band 
as a function of time during 2006-08. Each point represents one 
observation. Symbols as indicated. 
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between the two energy bands is that in the outburst 
peak, while in the soft band all components vary over 
similar strength (5-18%), in the hard band the 1 Hz 
component is clearly stronger (20-25%) than the other 
components. 

Figures 7 and 8 show the evolution for 2006-08 in 
the hard and soft band, respectively. The 0.1 Hz com- 
ponent is present, at similar frequency and amplitude in 
both bands, but there is only a single detection of the 
QPO (in the hard band), and the 1 Hz component is 
not detected. The 2006, 2007 and 2008 average power 
spectra of both bands are shown in Figure 9. In 2010 
there is no significant variability in either band. This is 
consistent with the source being in a softer, less variable 
state. 



Fig. 7. — Intensity, frequency and fractional rms amplitude of 
the components detected in the power spectra in the hard band 
as a function of time during 2006-08. Each point represents one 
observation. Symbols as indicated. 

The RXTE PCA has a larger effective area than 
the Swift XRT, so not all components observable with 
RXTE might be detected in individual Swift observa- 
tions. To investigate this, we averaged the power spectra 
over MJD 53552-53580 (hereafter the "outburst peak") 
in the 0.5-10 keV band, for a total exposure of ~14 ksec. 
One additional component was detected, at 1.5 Hz, but 
its single-trial significance is only 2.5o\ This component 
is similar to one that can be seen in the power spectra 
plotted in Z07 and RS07. We conclude that our Swift 
results in the 2-10 keV band are consistent with RXTE, 
although at reduced sensitivity, exactly as expected. 

We now look at the evolution of the power spectral 
components in 2005 in the soft band and compare it 
with the hard band. A representative power spectrum 
is shown in Figure 4, left panel. Figure 6 shows the 
light curve, evolution of frequency and the rms of these 
components with time. The most obvious difference 



Earlier studies (e.g., Homan et al. 2001, in the 
BHB XTE J1550-564) have shown that the QPOs are 
stronger and broader at higher energies. We investigate 
if the different components in the power spectrum have 
the same energy dependence based on the outburst-peak 
data. First we calculated power spectra separately in the 
hard and soft band and fitted these in two ways: with 
the v m ax and Q of the respective components "tied" 
(constrained to be the same) and only their integrated 
powers free, and with all the parameters free. The fits 
have a reduced x 2 of ~1.03 (with 11560 and 11554 
degrees of freedom, respectively), and similar values 
within errors for tied and untied parameters, except for 
the Q of the 1 Hz component, which in the tied fit is 
0.80±0.08, and in the untied fit 1.19±0.15 in the soft, 
and 0.66 ±0.1 2 in the hard band. The rms values of the 
QPO and the 1 Hz component are higher in the hard 
band than the soft band (both by ~9 %) in both tied 
and untied fits. This shows that the high frequency 
components are stronger at higher energies. The harder 
1 Hz component and QPO cause the extra power >0.4 
Hz reported above in Figure 3. The 1 Hz component 
is also broader in the hard band and perhaps should 
be fitted independently. Hence, we fitted the power 
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Fig. 9. — The rms normalized average power spectra for the years 
(top to bottom) 2006, 2007 and 2008 in hard and soft bands as 
indicated. The best fit model using one Lorentzian is shown in 
each frame. The power at high frequencies seen in some of the 
spectra above is not significant. 
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Fig. 10. — The correlation of QPO frequency with intensity dur- 
ing 2005-08. The energy bands are indicated. 

3.3. Correlations 

We finally look for correlations between power spec- 
tral parameters and intensity. We also looked at the 
correlations with hardness and, hardness being inversely 
correlated to intensity (see Figure 1), these are similar 
(but of course inverted relative) to the correlations 
with intensity. Hence, we show the correlations with 
intensity only. We first inspect the dependence of the 
characteristic frequency of the components on intensity. 
The 0.1 Hz component shows no correlation with 
intensity in either energy band at any time. As seen in 
Figures. 5, 6, 7 and 8 (stars) the component is present 
during 2005-08 over a wide range of count rates, but 
always at similar frequencies. The QPO frequency as 
a function of intensity is plotted in Figure 10. In the 
hard band (circles) there is some evidence of a positive 
correlation, as reported earlier by Z07 and RS07. In the 
soft band (crosses) there are few detections, and the 
dependence of frequency on intensity is not clear. The 
1 Hz component also has very few detections, and no 
clear correlation is seen. For each component, Q covered 
a wide range of values, as reported in Section 3.2. We 
inspected the dependence of Q on frequency, intensity 
and rms, but there were no clear correlations. 



spectra of all individual observations (discussed above) 
independently in both bands with all parameters free. 

Finally, to study the energy dependencies <2 keV, 
we calculated the outburst peak power spectra in the 
0.5-1 keV, 1-1.5 keV and 1.5-2 keV bands. All the 
outburst peak power spectra are collected in Figure 12 
and their best-fit parameters in Table 1. The 0.1 Hz 
component is detected with similar strength (11-15% 
rms), frequency and Q in all bands. The QPO is not 
constrained in the 0.5-1 keV band, it is ~8% rms in both 
the 1—1.5 keV and 1.5-2 keV bands, and nearly twice 
that amplitude again in the 2-10 keV band, indicating 
that this component is hard. The 1 Hz component is 
not constrained in the 0.5-1 keV band, is detected at 
~9% in the 1-1.5 keV band and more than doubles 
in amplitude to ^20% rms in the 1.5-2 keV and 2-10 
keV bands. So, the low frequency component is strong 
in all bands, while the high frequency components get 
stronger as we go to higher energies. 



The rms of the individual components for 2005- 
08 is shown in Figure 11 as a function of count rate. For 
both the 0.1 Hz component (stars) and the QPO (dia- 
monds) and in both energy bands, rms is anti-correlated 
to intensity. As there are no detections of the 1 Hz 
component (circles) at low intensities, its rms-intensity 
correlation is not clear. The 1 Hz component stands out 
in the hard band as the strongest of all components in 
the peak. 

4. DISCUSSION 

As this is the first time that power spectra are studied 
with XRT in such detail, we examined the behavior 
of the Poisson level and find it to be affected by 
instrumental effects (see the Appendix for details). We 
carefully compared our results with the earlier RXTE 
results using the 2-10 keV band, and found them to be 
consistent. We observe the Fast Rise Exponential Decay 
type light curve also reported in earlier works (C07, 
RS07, Z07 and ChlO). The HID does not resemble 
in either Swift or RXTE (Soleri et al. 2012) data the 
typical "q" -shaped track followed by many other BHBs. 
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Fig. 11. — The fractional rms amplitudes of the power spectral 
components as a function of count rate during 2005-08. The energy 
bands are indicated. Symbols as indicated. 



0.5-1 keV 




1-1.5 keV 




0.1 1 10 100 

■ 'Ccccy ( /) 



0.1 1 10 100 

Frequency (Hz) 



1.5-2 keV 




2-10 keV 




0.1 1 10 100 

Frequency (Hz) 



0.1 1 10 100 

Frequency (Hz) 

0.5-2 keV : 




0.1 1 10 100 

Frequency (Hz) 

Fig. 12. — The rms normalized average power spectra of the out- 
burst peak (MJD 53552-53580) in energy bands as indicated. The 
best fit model using three Lorentzians are shown. Corresponding 
parameters are presented in Table 1. 



The source remained in hard states throughout the 
outburst (C07, RS07, Z07, ChlO), until a transition to a 
softer state occurred in 2009 (Soleri et al. 2012), leading 
to our softest observations in 2010. 

The "type-C" QPO, the break frequency, Vb and 



Energy 
(keV) 


Vmax 

Hz 


Q 


rms (%) 


rms (%) 
0.03-10 Hz 


0.5-1 


0.23 ± 0.07 


0.08 ± 0.20 


11.9 ± 1.4 


13.0 ± 1.1 


1-1.5 


0.18 ± 0.02 
0.67 ± 0.05 
1.56 ± 0.10 


0.30 ± 0.10 
1.34 ± 0.79 
1.60 ± 0.72 


15.1 ± 1.0 
8.8 ± 2.1 
8.7 ± 1.3 


19.3 ± 1.3 


1.5-2 


0.12 ± 0.03 
0.72 ± 0.02 
1.07 ± 0.34 


0.45 ± 0.32 
3.91 ± 1.83 
0.0 (fixed) 


10.7 ± 2.2 
8.9 ± 2.3 

22.8 ± 1.5 


27.1 ± 1.3 


2-10 


0.18 ± 0.03 
0.69 ± 0.01 
1.79 ±0.23 


0.23 ± 0.15 
1.97 ± 0.36 
0.45 ± 0.20 


13.0 ± 1.3 

15.1 ± 1.5 
19.1 ± 2.4 


27.2 ± 1.1 


0.5-2 


0.21 ± 0.02 
0.69 ± 0.03 
1.6 ± 0.6 


0.08 ± 0.07 
1.80 ± 0.70 
0.91 ± 0.25 


14.3 ± 0.5 
6.3 ± 1.0 
9.5 ± 0.8 


18.7 ± 1.0 



TABLE 1 

Best fit parameters of the averaged power spectrum 
shown in Figure 12. all parameters were free. 



the hump frequency, v^, reported by Soleri et al. 
correspond to the QPO, the 0.1 Hz and 1 Hz compo- 
nents, respectively, in our analyses. We report these 
components, for the first time, in soft X-rays together 
with the harder band. The WK correlation between Vh 
and (Wijnands & van der Klis 1999) was reported 
by Soleri et al. between these components in SWIFT 
J1753. 5-0127. As there are very few simultaneous 
detections of these components in our data, we do 
not see a clear WK correlation in either band. The 
fractional amplitudes of all three components vary over 
similar range (5-18 %) in both bands, except the 1 Hz 
component in the hard band which is the strongest 
(20-25 %) amongst all components in the peak of the 
outburst. In the decay, all the components get stronger 
in both energy bands, with no detections of the 1 Hz 
component. The analysis in soft sub-bands and hard 
band (see Figure 12 and Table 1) in the peak of the 
outburst shows that the 0.1 Hz component displays a 
fiat rms spectrum i.e. the fractional amplitudes do not 
vary with energy. The QPO and the 1 Hz component 
cannot be constrained in the 0.5-1 keV band. They 
get stronger at higher energies, with the QPO getting 
stronger above 2 keV while the 1 Hz component gets 
stronger above 1.5 keV. 

The main new finding of our work (Figure 3) is 
that in the peak of the outburst (>110 c/s) the variabil- 
ity is weaker in the soft band than in the hard band; 
the difference is most pronounced at high (>0.4 Hz) 
frequency but also seen at low (<0.4 Hz) frequency. 
Previous work with XMM-Newton (Wilkinson & Uttley 
2009) had shown that in the decay, in 2006, the opposite 
was true: at low frequency the soft band was more 
variable than the hard band, while at higher frequencies 
both bands had similar fractional variability amplitudes. 
Our 2005-2008 results in the decay are consistent with 
this. Wilkinson & Uttley (2009) interpreted this extra 
soft-band variability in terms of a stronger contribution 
of the soft disk emission to the variability at low 
frequencies, which could be associated with accretion 
fluctuations intrinsic to the disk at these frequencies. 
Combining our new results on the outburst peak with 
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these earlier findings on the decay, we can say that 
in both frequency ranges the fractional rms spectrum 
becomes softer as the source decays, while instead the 
overall source spectrum becomes harder. 

These findings can be understood in terms of a 
simple two-component representation (hot inner flow 
and cool disk) of the picture of Lyubarskii (1997) and 
Churazov et al. (2001). The softer overall spectrum in 
the peak of the outburst is due to a relatively stronger 
disk contribution. That the fractional rms spectrum 
gets harder at all frequencies when the disk contribution 
increases can be explained if the observed disk variabil- 
ity has a lower fractional rms than the hot flow: the 
increased disk flux dilutes the variability at low energy. 
This would imply that while in the decay the disk at 
low frequencies is more variable (in fractional rms) than 
the hot flow, in the outburst peak it is less variable. So, 
we may be seeing in this outburst peak a transition to 
the less variable disk state that must also characterize 
the soft states. That in our outburst peak data the 
difference between low and high energy variability is less 
at low frequency can be due to a residual contribution 
of intrinsic disk variability at low frequency similar to 
that proposed by Wilkinson & Uttley (2009) and Uttley 
et al. (2011). 



It would be interesting to observe the variability 
behavior of this source when (if) it makes a transition 
back to the hard states, to see if it exhibits behavior 
similar to other BHBs during this transition, even if 
up to now it did not follow the usual "q" path in the 
HID. With our results, we demonstrate that the power 
spectral studies can be successfully performed with Swift 
XRT. The soft band provides very useful insight into the 
energy dependence of the variability which is essential 
to understand the contribution of the accretion disk. 
To understand better the role of the disk, variability 
studies of more BHBs with Swift XRT are indicated. 
We recommend caution while dealing with the Poisson 
spectrum of the power spectra generated with Swift 
XRT, however, as it has a different behavior than the 
RXTE Poisson spectrum and should not be assumed to 
be the same. 

ACKNOWLEDGEMENTS 

We would like to thank Kim Page for help in Swift 
XRT data reduction issues. This research has made use 
of data obtained from the High Energy Astrophysics Sci- 
ence Archive Research Center (HEASARC), provided by 
NASA's Goddard Space Flight Center, and also made 
use of NASA's Astrophysics Data System. 



APPENDIX 

The XRT is a grazing incidence Woltcr 1 telescope which focuses the incoming 0.3-10 keV X-rays onto a CCD 
of 600x602 pixels (23.6 x23.68 sq.arc-min). For a good time resolution (1.766 ms) and to mitigate pile-up in 
bright sources, the CCD is operated in the WT mode. In this mode, in wt2 configuration, only the central 200 
columns of the CCD are read out. Every 1.766 ms, 10 rows are compressed into a single row and then read out as 
a 1 dimensional image of 200 pixels. The spatial information in the other dimension is lost due to continuous clocking 2 . 

The arrival of photons follows a Poisson distribution. In a Leahy normalized power spectrum, we expect the 
Poisson noise level at 2.0. This holds true when the events are independent. However, in the case of pile-up, this does 
not hold true any more. Pile-up occurs when multiple X-ray photons incident on a 3 x 3 pixel region are read out as 
a single event of higher energy (Romano et al. 2006). The effects of pile- up on the Poisson noise spectrum have not 
been systematically studied yet for Swift XRT CCD in the WT mode. We examined the effects on the Poisson noise 
level of pile-up, which for the XRT CCD is expected to appreciably affect the data above 150 counts/sec. 

Power spectra were obtained using the same method as discussed in Section 2. We used observations of SWIFT 
J1753. 5-0127 covering count rates up to 190 c/s. Over the frequency range 50-100 Hz, where no variability is observed 
with Swift, we fit a constant to measure the Poisson level. Figure 14 shows the Poisson level as a function of count 
rate without, and for count rates above 150 c/s also with pile-up correction. It can be clearly seen that Poisson level 
decreases as the count rate increases. For the observations with count rates above 150 c/s, when no pile-up correction 
is applied, the Poisson level goes down to 1.86. If we apply pile-up correction to those observations, the Poisson level 
improves. However, it is worth noting that the Poisson level starts falling much below the nominal pile-up threshold 
of 150 c/s. So, even at much lower count rates, fixing the Poisson level at 2.0 is not good practice. This should be 
taken into account when fitting Swift XRT power spectra. We also investigated if there was any dependence of the 
Poisson level on photon energy, but no and obvious dependence was found. 

While inspecting the power spectra we also observed a power drop-off at frequencies above 100 Hz. Figure 13 shows 
the power spectra of observations of the sources SWIFT J1753. 5-0127 and MAXI J1659. 5-152 (another BHB) over the 
full frequency range and at higher frequency in the 0.5-10 keV band. The constant line at value 2 shows the expected 
Poisson level in an ideal case. A significant power drop-off is seen above 100 Hz, most clearly in the case of MAXI 
J1659. 5-152. This has been attributed to events that happen to fall across the 10-row boundary and are split into 2 
events of lower energy 3 . Again, this behavior is also seen in observations with count rates much below the nominal 
pile-up threshold. As it is not possible to identify and eliminate these events, in our analysis we put a cut-off at 100 Hz 



2 http://swift.gsfc.nasa.gov/docs/swift/analysis/xrt_digest.html 

3 see http : // 'www .swift.ac.uk / 'analysis / 'xrt / ' digest_cal .php#pow 
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Fig. 13. — Leahy normalized power spectra of two sources (as indicated) in the 0.5-10 keV band over the full frequency range and at 
higher frequency. Horizontal lines at power level 2.0 indicate the expected Poisson level in the absence of instrumental effects. 
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Fig. 14. — The Poisson noise level measured between 50—100 Hz in Leahy-normalized power spectra in the 0.5-10 keV band as a function 
of total count rate (uncorrected for pilc-up). The circles are observations that are not corrected for pile- up, while the filled squares are 
same observations above 150 counts/sec after pile-up correction. 



for all power spectra. 
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